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Abstract 

This paper presents a theoretical model for studying the dynamics of or- 
dering in alloys which exhibit modulated phases. The model is different 
from the standard time-dependent Ginzburg-Landau description of the evolu- 
tion of a non-conserved order parameter and resembles the Swift-Hohenberg 
model. The early-stage growth kinetics is analyzed and compared to the 
Cahn-Hilliard theory of continuous ordering. The effects of non-linearities on 
the growth kinetics are discussed qualitatively and it is shown that the pres- 
ence of an underlying elastic lattice introduces qualitatively new effects. A 
lattice Hamiltonian capable of describing these effects and suitable for carry- 
ing out simulations of the growth kinetics is also constructed. 
64.60.Cn,5.50.+q,82.20.Mj,64.60.My 
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I. INTRODUCTION 



One of the most fascinating problems in statistical mechanics is that of growth kinetics 
[IJ; the evolution of order following a quench from an initial high temperature state to one 
below the transition temperature. Theoretical studies of kinetics have mainly concentrated 
on understanding the behavior of models such as the time- dependent Ginzburg-Landau 
model or the kinetic Ising model || . Binary alloys have been viewed as convenient model 
systems for experimental investigation of these theoretical predictions KJ. Many alloys, 
however, exhibit complex ordered phases including long-period, modulated structures || 
which cannot be described by simple antiferromagnetic Ising models 0. The kinetics of 
ordering in these alloys should be significantly different from the usual models of ordering. 
These alloys can in turn provide models for both theoretical and experimental studies of a 
different class of kinetic phenomena. 

This paper presents model partial-differential equations for describing growth kinetics in 
alloys with modulated phases and analyzes the early stage kinetics. The kinetics of growth 
in the strongly non-linear regime have to be approached through solutions of the partial 
differential equations or through Monte Carlo simulation of a discrete lattice Hamiltonian. 
Such a discrete Hamiltonian, derived from the same atomistic description that leads to the 
partial differential equations, is also presented, and the qualitative aspects of this model are 
discussed. 

In recent years, there has been a concerted effort to understand phase stability in metallic 
alloys based on a realistic description of the "electron glue" responsible for metallic cohesion. 
One such approach, which has recently been applied to equilibrium phase transitions in the 
Cu-Au alloys with success |§||, is based on the effective medium theory (EMT) of cohesion 
in metals |H|. EMT provides a semi-empirical scheme for writing down the energy of a 
system of electrons and ions with an arbitrary configuration. This is its major advantage 
over first-principles approaches which have to rely on specific symmetries of the system such 
as a) completely random |L1] or b) ordered periodic structures |I2|| . By its very nature, 
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EMT is more approximate than the first-principles approaches; however, it does provide a 
model Hamiltonian which takes into account the metallic character of these alloys, describes 
the coupling of positional and configuration degrees of freedom |||9[] and can deal with size 
effects arising from the fact that the Ising spins are actually atoms which have different 



"sizes". 



The EMT expression for the total energy for a given configuration of ions provides the 
classical Hamiltonian describing the alloys |8|J§. This Hamiltonian is a function of both 
the positions of the atoms and their chemical identity. The latter can be represented by 
an Ising variable So far, our studies of alloys have been based on an EMT Hamilto- 
nian which includes the Ising variables and uniform lattice distortions, ie volume and shape 
distortions have been included but not volume conserving phonons From studies of 

compressible Ising systems, it is known that the coupling to uniform lattice distortions can 
change the nature of the Ising phase transition |14|,|15j. The most striking demonstration 
of this is an Ising antiferromagnet on a deformable triangular lattice where the system was 
found to order, through a first-order phase transition, into a striped phase accompanied by 
a lattice distortion which has the antiferromagnetic bonds shortened and the ferromagnetic 



bonds lengthened |nj. On a rigid lattice there is no ordering. Our studies of the ordering 
in CuAu, based on EMT with uniform lattice distortions, showed that including the lattice 
distortions leads to a first-order transition from the disordered, undistorted lattice to a "lay- 
ered" phase where the face-centered-cubic (fee) lattice has alternating Cu (spin up) and Au 
(spin down) planes stacked along the (100) direction and there is a tetragonal distortion (the 
Ll structure observed experimentally [[17]]) which shortens the antiferromagnetic bonds and 
stretches the ferromagnetic bonds ||||. This observation implies that the CuAu transition 
is in the same class as the triangular lattice Ising antiferromagnet coupled to uniform lattice 
distortions. Part of the physics entering the EMT Hamiltonian for Cu-Au can therefore be 
understood simply in terms of Ising antiferromagnets on frustrated but deformable lattices. 
The advantage of EMT is that the model is derived from a microscopic description of a given 
alloy and all the alloy-specific parameters are known. However, the mapping of the physics 



onto a simple model helps in understanding the physics of ordering in these alloys. 

The other ingredient in EMT which is different from the usual Ising model description 
of alloys is the size of an atom. This size is defined in EMT as that radius of a sphere within 
which each atom is neutral. The size is therefore a configuration dependent quantity which 
depends on screening in a metallic environment. It has been argued || that the modulated 
CuAuII phase [|17] is a consequence of such size effects and it will be shown here that these 
size effects lead to Ising models with competing interactions similar to ANNNI models 
The CuAuII phase has a one- dimensional modulation of the order parameter along a cubic 
direction and the period of the modulation, which varies as a function of composition, is ten 
lattice constants at the 50-50 composition. There is an alternate mechanism which could 
lead to modulated phases in metallic alloys, and this is driven by Fermi-surface instabilities 



ljj . EMT, in its present form, cannot describe subtle Fermi surface effects and may fail 
to describe modulated phases in alloys where this is the dominant effect. The model of 
kinetics that will be presented in this paper can encompass both classes of modulated-phase 
alloys. The model will be derived for size-effect alloys such as CuAu and the differences and 
similarities with Fermi-surface alloys will be discussed. 

II. MODEL OF KINETICS 



Pattern formation and kinetics of phase transitions are traditionally described by the 
time-dependent Ginzburg-Landau model 0. For a non-conserved scalar order parameter 
(?/>), this model Langevin equation can be written as, 

where J-"[ip] is the Ginzburg-Landau free energy functional, and f]{r,t) is a Gaussian noise 
term reflecting fluctuations of the heat bath [0,0: 

< r](r,t) > = 0, 

< r](r, t)7](r', t') > = k B T5(r - r')8(t - t') . (2) 



It should be noted that Fty] also acts as a Lyapunov functional; QS&£pM < 0. 

Most descriptions of kinetics in alloys have been based on the traditional ip 4 form of 
J-\ip\. In contrast, our discussion of kinetics will be based upon a Fty] derived from the 
microscopic EMT Hamiltonian. 

We use as our prototype alloy the CuAu (50-50) alloy which exhibits modulated phases 
and, in equilibrium, has a first-order phase transition accompanied by a lattice distortion. 
The latter shows the relevance of the coupling to phonons and uniform lattice distortions 
and raises questions about the description of metastable and unstable states and the concept 
of a spinoidal. The presence of two first order transitions, from the disordered to the 
modulated CuAuII phase, followed by the transition from CuAuII to the zero temperature 
unmodulated CuAul structure MM, leads to unusual quench-temperature dependence of 



the kinetics and a Swift-Hohenberg-like |18| description of the early-stage kinetics. 

The free-energy functional for the CuAu alloy has been constructed by numerically 
calculating the mean-field free energy from the EMT Hamiltonian §| and we will only briefly 
discuss the procedure and the results in this paper. The Hamiltonian includes uniform lattice 
distortions and the Ising degrees of freedom, H({e a }, {sj}). Where Sj are the Ising variable 
on lattice sites i, and the {e a } are the set of uniform lattice distortions. The partition 
function for this model can be written as 

Z = J U a de a J Xtyexp(-F[V,e a ]) (3) 

ie, one integrates over the lattice distortions and the order parameter field ip. The Ginzburg- 
Landau functional J-"[ip] is obtained by making a saddle-point approximation to the e a 
integrals, ie minimizing F with respect to these variables. In our numerical procedure, 
the mean-field free energy corresponding to a given set of ip and a given set of e a was 
calculated from the EMT hamiltonian and this function was then minimized with respect 
to the e a . This is equivalent to making the magnetothermomechanics [H| approximation. 



The minimization of with respect to ip leads to the mean-field description of the CuAu 
ordering transitions [f|. In this paper, this Ginzburg-Landau functional is used to describe 



the kinetics of phase transitions in these alloys. 



The form of for CuAu, is most conveniently written in momentum space [2(| 

T[^\ = £[a(T - T ) + u(q)W q ^ q 
Q 
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+ ^(E ^q 1 ^q 2 ...^qj(q 1 + ...q 6 )} • (4) 



where 



f (g) = e g 2 2 + (gi - go) 2 - ?o • 

Here q± is the magnitude of a wavevector in the plane perpendicular to the CuAu ordering 
which is taken to be in the ^-direction. 

Before discussing the model of kinetics that follows from this Ginzburg-Landau func- 
tional, we briefly discuss the applicability of this model to equilibrium phase transitions in 
CuAu. Because of the form of u;(q), the Ginzburg-Landau functional, .F^], is minimized by 
a one-dimensional periodic function, which describes a configuration comprised of a periodic 
array of antiphase boundaries separating regions with ip — ±1. This is the source of the 
CuAuII phase. The appearance of stationary one-dimensional, periodic patterns is also a 



feature of the Lyapunov functional associated with the Swift-Hohenberg equation [TJ) . The 
isotropic form of u(q) implies that the one-dimensional pattern is rotationally invariant in 
both models. However, the CuAuII structure observed in Monte Carlo simulations |||§, 
based on EMT, and in experiments, has the modulation locked along one of the cubic axes, 
and there are lattice distortions accompanying this transition. The unusual fourth order 
term (3 rd term in Eq. (||)) arises from the coupling to uniform lattice distortions. This term 
is responsible for the CuAul type ordering MM which, as discussed earlier, is in the same 



class as the Ising antiferromagnet on an elastic triangular lattice [|16||. This type of ordering 



involves two broken symmetries, one corresponding to the twofold Ising degeneracy and one 



corresponding to the three possible strain directions on the lattice. The CuAuII ordering 
involves a lattice distortion along the direction of the one-dimensional modulation, and an 
additional symmetry breaking in the choice of this direction. It is expected that the coupling 
to the lattice is the source of the commensurate CuAuII ordering. In the Landau theory, 
presented earlier || , the modulation was assumed to be along one of the cubic directions. 
The effects leading to the choice of this direction, within mean-field theory, are currently 
being investigated. It is clear from our Monte Carlo simulations that the EMT hamiltonian 
can describe the lattice distortions accompanying the modulation ||||. 

The function u(q) can be measured in diffuse scattering experiments [f2Tp2f . The EMT 
model would predict an isotropic distribution in the plane perpendicular to the ordering 
direction (c/form of ou(q) in Eq. (f|)) and this isotropic distribution would have a peak 
displaced from the CuAul superlattice Bragg peak by a distance go- Experimental mea- 
surements on CuAu show an essentially isotropic distribution with a peak at go — 27r/10a 
where a is the lattice constant plfl . The EMT estimate for g is ~ 2ir/Aa ||. This was, 
however, a very crude numerical estimate which was only meant to show that go is nonzero. 
The essentially isotropic distribution of CuAu is in sharp contrast with diffuse scattering 
intensities in Cu^Pd [S2| where the distribution has well defined peaks along the q x and q y 



directions. These peaks have been explained on the basis of a Fermi-surface mechanism p3 
The essentially isotropic distribution in CuAu suggests that a different mechanism is respon- 
sible for the modulated phase. The small anisotropy seen in experimental diffuse scattering 
measurements in CuAu would be difficult to obtain from our numerical calculations based 
on EMT. It could also be that there is a Fermi surface contribution which is missing from 
EMT that makes the distribution anisotropic. Because of the small anisotropy, it is reason- 
able to investigate the isotropic model and discuss the consequences of the anisotropy. The 
isotropic model is interesting in its own right since it implies an the breaking of a continuous 
symmetry at the CuAuII transition, where the system chooses a direction of modulation. 

The above discussion shows that J 7 ^] is a good model for describing kinetics of these 
alloys. The kinetic equation derived from the model ^ r [V ; ], is 



^ (r ' t] - 1 [ a(T - T ) - qt] - ^ + {ql + Vi) 2 } V»(r, t) - ^ + ^(r, f), (5) 



where 

III. EARLY STAGE KINETICS 

The early stage kinetics of CuAu, following a quench, is defined by a linear stability 
analysis of the disordered, if) = state. This implies neglecting the 8f/8ip term in 
Eq. (Bp. The linear theory is very different from the usual model of a non-conserved 
order parameter The difference is evident from an analysis of the solution to the 



linearized equations without the noise term. The easiest route to this solution is through 
transformation to momentum space where the different modes decouple and the solutions 
for if>(q) are, 

if;(q,t) = exp(D(q)t) (6) 

where D(q) = — a(T — T ) — uj(q). This solution shows that modes with D(q) greater 
than zero grow exponentially with time. Eq. (|6|) is the generalization of the Cahn-Hilliard 
solutions describing continuous ordering p5| to alloys with modulated phases. If none of 
the eigenvalues, D(q), are greater than zero then there is no continuous ordering and the 
disordered state is either stable or metastable. In contrast to the Cahn-Hilliard solutions, 
the fastest growing mode has a finite wavevector, q Q and the growth rate depends only on 
the magnitude of the wavevector and not its direction. This is reminiscent of Cahn-Hilliard 
description of spinodal decomposition but there the wavevector of maximum growth depends 
on the temperature |T],|25| whereas in our model, the linear dynamics chooses a unique length 
scale which is a characteristic of the system. Since the interesting pattern formations take 
place in the plane perpendicular to q z , we will analyze the equations in the q z = plane. This 
linearized CuAu equation is then identical to the linearized version of the two-dimensional, 



stochastic Swift-Hohenberg equation [26 



di/j(r, t) 



{ e -(l + Vi) 2 }^(r,t), (7) 



dt 

if the control parameter e is identified with go 4 — a(T— To) and the natural wavevector go is set 
to 1. The Swift-Hohenberg equation is used to model pattern-formation in Rayleigh-Benard 
convection and has been studied extensively both in one and two dimensions [ f26|j . 

The linear dispersion relation (D(q)) is plotted in Fig. (1) for various temperatures. 
Three different temperature regimes can be identified. For T > Tq + qo 4 /a, all modes de- 
cay and there is no continuous ordering, this is the regime where the disordered state is 
metastable and ordering takes place by nucleation. For (T + go 4 ) < T < T , a band of 
wavevectors centered around go grows and at T = To, the q = mode becomes unstable. 
The q = structure corresponds to the unmodulated CuAul phase. In the Cahn-Hilliard 
description of spinodal decomposition, the growth rate at q = is always zero, and in its 
description of continuous ordering, the q = mode always has the fastest growth rate. It 
should be remembered that the linear dispersion relation is isotropic in the plane perpen- 
dicular to q z and all modes with a given magnitude of q grow at the same rate. This is 
expected to give rise to an interesting morphology of the early stage ordering process [E7[ . 



Experimental sudies of kinetics are based on the structure factor. The structure fac- 
tor obtained from the linearized equation at a temperature within the second temperature 
regime is shown in Fig. (2a). The plot shows the evolution of the structure factor as a 
function of time. The peak of the structure factor remains stationary CIS CIS db small band 
of wavevectors centered around this peak grows as a function of time. The modes outside 
this band are seen to decay with the growth of order. At a lower temperature, the structure 
factor at q = would start to grow, but the growth rate predicted by the linear theory 
always has its maximum at q = go- The inclusion of the Gaussian noise term changes the 
form of the structure factor (generalization of the Cahn-Hilliard-Cook theory p5| ) and can 
be written as 

S(q, t) = S (q) exp(2D(q)t) + ^L.(exp(2D(q)t) - 1) (8) 
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where So(q) is the structure factor in equilibrium at the temperature prior to the quench. 
The inclusion of the Cook term introduces minor modifications of the structure factor which 
are illustrated in Fig. (2b). These plots show that an experimental investigation of the 
early stage kinetics of CuAu in this temperature regime should show a peak growing at a 
finite wavevector (measured with respect to the superlattice Bragg peak) and there should 
be a ring of wavevectors at which the structure factor is a maximum. The question of 
whether the early-stage description is ever valid in these systems can only be answered 
after a detailed study of the complete non-linear equations or from Monte Carlo simulations 
of the appropriate model Hamiltonian. It has been argued that the length of time over 
which the early stage kinetics remains valid is given approximately by 1/D(q ) and depends 
logarithmically on the range of interaction f28f . In our model, we have an effective infinite- 
range four-spin interaction arising from the coupling to the lattice distortions and this might 
imply that the range of validity of early-stage kinetics will be significantly enhanced in these 
systems. The evolution of the lattice distortion as a function of time is a novel feature of 
these systems and experimental investigations should shed light on the coupling between 
the evolution of order and the evolution of the lattice distortions. 



IV. DISCUSSION OF NONLINEAR EFFECTS 

The effect of nonlinearities may change the kinetics in various ways. In the simplest 
case, the nonlinearities mainly come in to dampen the exponential growth. Some simple 
considerations of the nonlinearities in our model indicate that there can be more significant 
differences between the linear and nonlinear dynamics. The periodic solutions to which the 
■0 = state becomes linearly unstable have wavevectors lying in the range [ki, fo], where 

fc 1>2 = (g 2 ± vV-^T-To))) 1 ^ (9) 

For temperatures close to but less than T\ = To + qo 4 /a, the width of the band of wavevectors 
grows as |T — Tip 2 . As the temperature approaches T , the q = mode becomes unstable 
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and below To, the width grows as (T — To) for T close to To. This band of wavevectors, chosen 
by the linear dynamics, usually define the range in which the stationary periodic solution 
is found at a given temperature; ie, one could examine the periodic solutions lying in the 



range [fo, ki], and the solution with the lowest free energy becomes the selected pattern p9 
However, because of the negative fourth order term in •T[V ; ], there is a first order transition 
from the disordered to the ordered periodic phase. This implies that, in a given temperature 
range, periodic solutions which are outside the band predicted by the linear dynamics could 
arise through a nucleated process and could become the stationary solutions. The condition 
that a order parameter of the form ip(r) = ipexp(iq ■ r) is a solution to 

is 

m > -^2! 

compared to the condition for linear stability, D(q) > 0. For example, the q = state 
becomes the stable stationary solution at a temperature higher than its instability temper- 
ature Tq. The role of nonlinearities is different in this model than in the Swift-Hohenberg 
model. Since the negative fourth order term is a consequence of the coupling to the lattice 



fqhqj , the investigation of non-linear dynamics has to take proper account of this coupling. 
The other crucial difference between this model and the Swift-Hohenberg model, is that the 
underlying lattice in CuAu introduces a second length scale comparable to the length scale 
chosen by the linear dynamics (27r/g ) and may be instrumental in locking in the direction 
of modulation. These effects are , in principle, all there in the nonlinear terms of our model 
equations if we are careful in keeping the Umklapp terms involving nonzero reciprocal lattice 
vectors in Eq. (^). However, it might be more profitable to investigate these effects directly 
in an effective lattice model. 

Numerical solutions of the full non-linear Swift-Hohenberg equations have shown that 
the lamellar patterns exhibit phases which are akin to nematic and smectic ordering in liquid 
crystals |26]. In CuAu alloys, the lamellar patterns have a length scale comparable to the 
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underlying lattice and an analogy maybe drawn between these systems and nematic ordering 
in the presence of an underlying lattice JJIJ . The late stage growth laws in Swift-Hohenberg 
systems is expected to be different from usual models of non-conserved order parameters 



since the wavevector defining the modulation has a continuous symmetry |56fl . It would be 
interesting to study the late stage growth laws in our model which has the extra feature of 
coupling to the lattice. 

V. LATTICE HAMILTONIAN 

We have argued that the study of non-linear dynamics is CuAu-like alloys may be best 
approached through a discrete lattice Hamiltonian. With this in mind, a lattice Hamiltonian 
has been constructed which has the salient features of the EMT Hamiltonian when applied 



to ordering on a deformable lattice |31 |. The model is an Ising antiferromagnet on an elastic 
lattice with an unusual pair interaction term which arises from the size difference between 
the two chemical species making up the alloy. The source of this term in EMT is the density 
dependent cohesive energy [|K],(J]. When expanded in terms of the spin variables, a part of 
this energy gives rise to a term in the Hamiltonian which is an antiferromagnetic interaction 



between two atoms which share a common nearest neighbor |52 |. This interaction depends 
on the bond lengths and therefore can be different along different lattice directions in the 
presence of deformation. The strength of the antiferromagnetic interaction depends on the 
size-difference of the two atoms and is zero for atoms of the same size. This model can be 
written as, 

H = EA 1 ~ ee ")£<y >a s ^ + E ap M a ,e a e, + £,(£ <ii>a %) 2 • (10) 

The parameters can all be obtained from EMT, given a specific alloy, however, we will 
treat this as a general model with arbitrary parameters. This should then form a generic, 
minimal model for alloys where Fermi surface effects do not dominate. The first two terms 
in Eq. ( ]TD| ) describe an Ising model coupled to uniform strain through the coupling con- 
stant e and an elastic Hamiltonian defined by the tensor M a p. This model for J > has 
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been analyzed on the triangular lattice ]16[ and gives rise to a first order transition from a 
disordered to a striped phase. The last term in the Hamiltonian is the unusual term arising 
from size effects. Since K > 0, this term favors a spin configuration where, for each site, 
the sum of the nearest neighbor spins along a given direction add up to zero. In terms of 
atoms, this says that when atoms of different sizes are present, it is favorable to have a big 
atom and a small atom as nearest neighbors. Neither an antiferromagnetic, nor a striped 
phase meets the optimum demand of the K term and a mean-field argument shows that, 
on a triangular lattice, modulated phases occur beyond a critical value of K [fffil . This is 
reminiscent of the ANNNI model. The one crucial feature which distinguishes the current 
model from the ANNNI model is that all order-disorder transitions are first-order because 
of the coupling to the lattice and makes it a more attractive model for describing alloys. 
Qualitatively, the source of the modulation can be inferred from looking at the striped phase 
of the triangular lattice antiferromagnet which has alternating chains of up and down spins 
and the sum over nearest neighbor spins for each direction adds up to two. If a modulation 
is introduced along the ferromagnetic direction, such that after N spins, the spin up and 
spin down chains are interchanged, then the spins at the domain boundaries have nearest 
neighbor configurations which do satisfy the K term requirement. The reason for the choice 
of the ferromagnetic direction has its origin in the preferred lattice distortion. Therefore 
in spite of the K term being isotropic, a unidirectional modulation occurs. This scenario 
is very similar to the CuAuII ordering and investigations on the fee lattice are currently 
underway. At this stage, we would like to present this model as an interesting variation on 
competing interactions models and suggest that the study of kinetics of this model would 
lead to novel features which can be tested through experimental studies of alloys. 

VI. CONCLUSION 

As mentioned at the beginning of the paper, there is a class of modulated alloys where 
the modulation is due to Fermi-surface effects P3J. In these alloys, the q— dependence of the 
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second order term in is not isotropic; ie, u(q) has well defined minima along directions 
chosen by the Fermi-surface nesting vectors. The linear dynamics in these alloys chooses 
a length scale and a set of directions and the symmetry breaking that takes place is not 
continuous but discrete, Ising like. Also, the role of the lattice distortions are expected to be 
minimal and therefore the nonlinearities may play a simpler role and the late stage kinetics 
may follow the same scaling laws as in the usual models. 

A complete description of ordering kinetics in alloys has to take into account the effect 
of phonons and local lattice distortions. Allowing for fluctuations in bond lengths also 
leads to an additional size-effect derived term in the lattice Hamiltonian. Under certain 



circumstances, this can give rise to a long-range pair interaction |33|. The EMT prediction 
for this strain-mediated pair interaction term and its effects on kinetics is currently being 
investigated. 

In conclusion, we have presented a model for describing kinetics of ordering in alloys 
which have modulated phases. The linear theory shows that the kinetics is similar to that 
of Swift-Hohenberg models, but the nonlinear terms are distinctly different from the Swift- 
Hohenberg equations. A lattice Hamiltonian which captures the essential features of the 
microscopic model of these alloys has been constructed and simulations based on this model 
should be helpful in understanding both the statics and dynamics of ordering in these alloys. 
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FIGURES 

FIG. 1. The linear dispersion relation D(q) vs. \q±\. The dispersion is isotropic in the |qj_| 
plane. Four temperatures are shown with the bottom most curve at the temperature T\ = To+go 4 / 
and the topmost curve at the temperature To. 

FIG. 2. (a) The evolution in time of the structure factor obtained from the linear theory without 
the noise term. The plots show S(q) at a temperature T such that T\ < T < Tq. The crosses are 
at the earliest time and the squares at the latest time, (b) The evolution in time of the structure 
factor predicted by the linear theory with the noise term included. The parameters used are those 
appropriate for CuAu. 
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